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Abstract 

In this paper we study stochastic dynamics which leaves quantum gravity equilibrium distribu- 
tion invariant. We start theoretical study of this dynamics (earlier it was only used for Monte-Carlo 
simulation). Main new results concern the existence and properties of local correlation functions 
in the thermodynamic limit. The study of dynamics constitutes a third part of the series of papers 
where more general class of processes were studied (but it is self-contained), those processes have 
some universal significance in probability and they cover most concrete processes, also they have 
many examples in computer science and biology. At the same time the paper can serve an intro- 
duction to quantum gravity for a probabilist: we give a rigorous exposition of quantum gravity in 
the planar pure gravity case. Mostly we use combinatorial techniques, instead of more popular in 
physics random matrix models, the central point is the famous a = — | exponent. 



'Postal address:INRIA - Domaine de Voluceau, Rocquencourt, BP105 - 78153 - Le Chesnay Cedex, France. 



1 



Contents 



1 Introduction 



2 One dimensional gravity 

2.1 Equilibrium distribution . . 

2.2 Gravity as a queucing mode] 



3 Spin Complexes 

3.1 Cell structures! 




3.1.1 Abstract complexes 

3.1.2 Topological complexe 

3.1.3 Cell surfaces .... 






3.1.4 Local observer (root) 

3.1.5 Mnves 






3.1.6 Automorphism group 


3.2 


Metrics and Curvature 



4 


Equilibrium planar pure gravity 


4.1 Definitions and combinatorial approach 




4.1.1 Uniform asymptotics 




4.2 RMT approach 



5 Linear boundary dynamics 

5.1 Local Pure Growth 

5.2 Reversible boundary processes 



6 Nonlinear boundary dynamics 

|6.1 Quadratic quasi processes 

|6.2 Generating functions 

|6.3 Analysis of the functional equation 

|6.4 Correlation functions 

|6.4.1 Tree representation 

6.4.2 Local curvature 



7 Internal dynamics 



1 Introduction 



Some history I should say some words about the history of discrete gravity. Classical gravity deals 
with a smooth (not necessarily four-dimensional) manifold M, pseudo- metric tensor gij on it and the 
classical Einstein-Hilbert action 



S= (XR(x) + L(x) + / u)y / det gdx 
Jm 



the various stationary points of which are studied. Here R(x) is the intrinsic curvature at the point 
x, L - some functional of matter fields <j>(x). In the pure (no matter) gravity case L = 0, we consider 
only this case here. Quantum gravity takes into account not only stationary points but also all other 
configurations with some weights, that is with a formal (but which becomes positive for Euclidean 
metrics) density 

Z~ l exp(—fiS) 

on some configuration space 0, of matter fields and metric tensors. All earlier attempts to do this 
brought the conclusion that f2 should also include smooth structures on M and even M itself, that 
is the space, its topology, should be random. Now the only reasonable way to pursue this program 
is to discretize everything from the beginning and then to perform some scaling limits. That is the 
space becomes a finite complex, smooth structure becomes a piecewise linear structure, metrics and 
curvature are encrypted in one-dimensional and two-dimensional skeletons of the complex, matter fields 
are spins which live on the cells of the complex. It appears that such quantization (discretization) 
is equally applicable to other physical systems: relativistic particles, strings etc., but with different 
interpretations. For example, the quantized (in such a way) string consists of a two-dimensional 
complex (representing a coordinate system and metrics on the string itself) and spins - vectors in R d 
which provide a mapping of the vertices of the complex into the (i-dimensional Euclidean space, thus 
approximating the classical string. 

The discretization of the classical gravity was first considered by Regge where he gave defi- 
nitions of some exact mathematical objects related to the classical general relativity: finite discrete 
space time, its curvature and Einstein-Hilbert action. It was afterwards included in the fundamen- 
tal monograph |3(| but in seventies it was still considered outside of the main streamline of physics 
and only rare papers were devoted to it. Among them however there was a well known paper by S. 



Hawking [ 32 1 where the applications to quantum gravity were discussed. 

In eighties there are already more than 100 papers concerning discrete quantum gravity. 

In nineties the the number of papers is more than 1000 and still grows at the moment. Mainly it is 
due to the appearing algebraic formal techniques to deal with such problems. This formal techniques 
follows physical insights on relations of quantum gravity with string theory random matrix models 
etc. Moreover, recent papers in theoretical physics often contain the following sententions: "Two- 
dimensional random geometry is now placed at the heart of many models of modern physics, from 
string theory and two-dimensional quantum gravity, attempting to describe fundamental interactions, 
to membranes and interface fluctuations in various problems of condensed matter physics" , see [^] . 

For a probabilist the quantum gravity is a source of inspiration and also new mathematics and new 
philosophy of probability. The paper can serve an introduction to quantum gravity for a probabilist: 
it is a mathematical text on the quantum gravity for the planar pure gravity case. 



Dynamics contre equilibrium Mostly we consider combinatorial techniques, instead of more pop- 
ular in physics random matrix models, the central point is the famous a = — | exponent. Another 
goal of the paper is to consider stochastic dynamics which leaves quantum gravity equilibrium distri- 
bution invariant. We start theoretical study of this dynamics (earlier it was only used for Monte-Carlo 
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simulation). The study of dynamics constitutes (but mainly it is self-contained) a third part of the 



series of papers (see [43, 44]) where more general class of processes was studied. These processes have 
also some universal character in probability: they cover most concrete processes. Also they have many 
examples in computer science and biology. 

Here the probability is the classical probability. The quantum gravity constitutes a bunch (a 
lot !) of papers overfilling last 10 years well-known physical journals. Discrete quantum gravity is 
now considered as a promising direction towards unifying largest and smallest scales in the nowadays 
picture of nature. I consider one part of this field which evidently uses probabilistic intuition but it 
is difficult to find even formulations (I do not mention proofs !) which could be satisfactory for a 
mathematician: even when the "probabilities" are hopefully positive they are not normalized. And 
this is not because of negligence of the authors but because some deep reasons seem to be behind the 
curtains. 

In the existing physical literature a permanently developing algebraic and geometric techniques 
overwhelms the subject. Thus it can be useful to step away from algebra and geometry, discussing 
some simple probabilistic aspects of quantum gravity: even such simple project appeared to rise many 
natural but still not answered questions. 

There are now two variants in the discrete approaches to quantum gravity: Quantum Regge Cal- 
culus (where links (edges) have lengths as random variables) and Dynamical Triangulations (where 
lengths of edges are constant). The word dynamical in the second approach is a little bit mislead- 
ing because there is no dynamics at all in this approach: main techniques uses Gibbs equilibrium 
distributions on large matrices. That is why I will call here these approaches equilibrium. 

The dynamics appeared earlier in Monte Carlo simulations of quantum gravity. Here I try to give 
a probabilistic (not numerical) study of relevant Markov processes. What is new here (I do not know 
earlier rigorous results) is that we want to advocate not numerical but analytic and probabilistic studies 
of such processes. Why such processes can be useful not only in computer Monte-Carlo experiments 
but also as giving theoretical information ? There are many reasons - we give here a short list. 

• Well known difficulty in averaging over all topologies is that, in 4 dimensions, it includes some 
questions which are known to be algorithmically unsolvable. Dynamics substitutes this problem 
with a new one: instead of averaging we are looking for a process (with arbitrary initial state) 
which will generate all topologies it can generate. This process should have some symmetries 
but also it should be a legitimate (for example non-exploding) stochastic process. 

• What one would like to have (as in the stochastic quantization in quantum field theory and 
Glauber dynamics in statistical mechanics) is a Markov process leaving Gibbs measure invariant. 
This is quite natural in quantum field theory where there are Whiteman axioms and in statistical 
mechanics where there is a deterministic dynamics more fundamental than the Gibbs measure 
itself. In quantum gravity both these factors are absent and an alternative viewpoint could be 
advertised: that the process itself can be taken to be more fundamental than the Gibbs measure 
itself. 

• Dynamics allows to consider the region below the critical point where equilibrium distribution 
has no sense. On the contrary this region is even more natural for the dynamics - like a growing 
universe (in the computer time, the term which I know from a paper by A. Migdal). Moreover 
dynamics gives also some sense to distributions in the critical point without performing scaling 
limits. I do not know the physical counterpart of all this but its naturalness from probability 
point of view is evident. 

• I have absolutely no physical arguments for the choice and even relevance of the dynamical 
models, but that is also true for all modern approaches due to the lack of experimental confir- 
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mation. The leading thread can only be probabilistic intuition and beauty. Relevant question 
are: what is universality and generic situation ? It was argued recently, see |pq1 , that computer 
science could play some role in future physical theories. Probabilistic aspects which we discuss 
here make this relation quite evident by a preliminary model of the universe growing via some 
grammar (more exactly a graph grammar) similarly to the random evolution of a language. 

Mathematical thermodynamic theory existing for statistical mechanics and quantum theory 
brought many new ideas. It is some surprise that an attempt to construct similar theory for 



growing complexes brought quite unexpected phenomena (see |43|] , J|4|] ) (hopefully having some 
physical significance). One of the effects is that one cannot fix an origin in an infinite universe 
without Zermelo axiom, any constructive introducing of a local observer changes drastically the 
space time in his neighborhood. 

It cannot be easy to find critical exponents by Monte-Carlo simulation because the asymptotic 
is dominated by the exponential term which depends strongly on the details of the model. What 
is usually simulated is the uniform distribution on the set of triangulations with fixed number 
of cells. If we consider a growing complex then we could not find a Markov process giving the 
necessary exponents (the famous a = — |) but only some random transformation of measures, 
that is called usually a nonlinear Markov processes, giving these exponents. 



Contents of the paper One-dimensional case (section 2) is useful in particular as emphasizing 
links between classical probability and two-dimensional quantum gravity. 

In section 3 the minimum of necessary definitions are given concerning complexes and curvature 
in two dimensional case. 

Section 4 contains introductory definitions, problems and some known results. In section 3.2 we 
give a short exposition of RMT approach to pure planar gravity, the only goal of this exposition is to 
emphasize some points, related to the combinatorial approach. 

In section 5 we study a dynamical model, where the cells are appended in random to the boundary 
of the disk. This model is solvable (via random walks) and we calculate some main quantities. The 
exponent for this model is —2 and thus it belongs to a different universality class than models accepted 
in physics. But continuum limit in this model is well defined, gives space with a constant curvature 
as in the physical model. 

Section 6 is the central in the paper. We construct and study nonlinear Markov processes (where 
also changes are possible only on the boundary) which render the equilibrium distribution invariant. 
We use the Tutte functional equation method to prove that one gets — | exponent. We develop new 
combinatorial techniques to study local correlation functions. 

In section 7 we consider dynamics where changes can be done elsewhere in the complex. We study 
large time behavior of such Markov processes. 



Acknowledgments I thank L. Pastur for elucidating to me some points of the Random Matrix 
Theory and S. Shlosman for reading the paper and very helpful comments. 



2 One dimensional gravity 

For the physical interpretation of the one-dimensional gravity and many beautiful calculations we refer 
to chapter 2 of Ambjorn's lectures Our goal here is to give a probabilistic viewpoint and discuss 
new approaches. There is no topology in one dimension: the underlying structure (cell complex) is 
one-dimensional - a linear graph. 
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Figure 1: Linear spin graph 

A chain of symbols from some alphabet can be considered as a function on vertices of such graphs 
(see Figure [j]) . 

2.1 Equilibrium distribution 

We give now the basic definition in more abstract terms than in |HJ , without prior embedding in Eu- 
clidean space: this corresponds more to Polyakov string quantization. We shall consider distributions 
on the set of finite linear spin graphs (sometimes we use the terminology from our previous papers 
but mostly it can be skipped). That is the distributions on the set O of strings (here string comes 
from computer science terminology) a = {s\, sjy}, where N = 0,1, 2, ... and Sj £ S, where S is an 
alphabet (spin space). For example S can be the unit sphere in Z d or in R d . Case N = corresponds 
to the empty string with no 5-value prescribed. Define the nonnegative measure on £1 by 

N-l 

Q(a) = exp(-/iiV - /3 ^(/x( Sm ) + f 2 ( Si , s i+1 )) (1) 

8=0 

for some function /. It is convenient to assume that some element sq is fixed. Simplest example is 
when s G Z d , fa = 0, fi(s) = oo for all s except a finite set. This measure can be normalized if 

a) < oo 

17 

Intuitively, a sequence of arrays sq, so + s\, sq -Mi + ■•• + sjv can be considered as a "random walk". 
But it is quite different from the classical random walks. We shall see below its relationship with some 
computer science problems. We shall also see how this formal object can be tied to the Euclidean 
space: in physical papers one can also see similar steps - abstract object (random triangulation, 
internal metrics etc. ) mapped finally to the physical space-time. 

Simplest examples In the first example S = Z d , fa = 0, fi(<r) = 1 if \a\ = 1 and oo otherwise, 
assume sq = 0. Otherwise speaking we have the non-normalized distribution exp(— l(s)) on all possible 
finite paths r = (0, s\, s± + S2, —, «i + «2 + ••• + sn) i n Z d , starting from where iV = l(r) is the length 
(number of steps) of r. It does not always exists. There exists < [i CT < oo such that the series 

Z = J^exp(-//Z(r)) 

s 

for the partition function converges for \x > fi cr and diverges for \x < fi cr . In our case fi cr = In 2d. 

For the second example S = R d ,e~^ 1 ^ = 5(\s\ — 1), fa(s,s') = 4>((s,s')) for some bounded 
function <fi of the angle between the two vectors. For the third example S = R d , fa = 0, fi(s) = s 2 = 
(s 1 ) 2 + ... + (s d ) 2 . 

These examples, highly simple and having nothing special from the probabilistic viewpoint, cor- 
respond to one-dimensional analogs of rather famous actions: free relativistic point particle in d- 
dimensional space-time (I - length parameter in the Euclidean space, L = L(0, x) is a path from to 
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x) 

S(L) = n j dl 

Hilbert-Einstein action (k is the curvature of the curve embedded in an Euclidean space) 



S(L) =n f dl + X f 
Jl Jl 



k\ dl 



' L JL 

and bosonic string action {g is some metrics on the parameter interval) 

72, 



s{L) = -,{! vsmoKfl^+M) 



One can consider the introduced distribution as a quantization of the corresponding classical action. 
Each link of the discrete path is assumed to have unit length and thus the length of a path is the 
number of links. 



We shall discuss only the first example, other two are similar, see [31]. For /i > /x cr we can define 



^ N .\( X ) e (-V + Vcr)N 



the probability distribution on the set of all finite paths starting from 

P( r ) = Z _1 exp(-Z(r)) 

Green function are defined as a measure on Z d 

G(x) = Yl exp(-/iZ(r)) = ^C(iV;x)exp(- / uiV) 
s.O— *x 

N 

where C(N; x) is the number of paths from to x of length N, P^y(x) - iV-step transitions probabil- 
ities from to a; for the classical simple random walk in Z d . The number of such paths (by the local 
limit theorem) is 

C(N;x) ~ (2vr)-f(2d) iV ' N~2 

if x = 0{^/N). Green functions have their origin in physics and they also look like Green functions 
for a Markov process 

where L are all paths from to x. But there is no Markov process here. 
The following observations are important: 

• neither fi cr nor the exponent — f do not depend on the choice of x. 

• [i cr is not universal, it depends on the dimension and on the lattice. Also we will get different 
values of n cr if we take piecewise linear paths in R d with sides of fixed length (see pll] ) : 



however the exponent — | does not depend on the lattice, which follows immediately from the 
local limit theorem; 

for = [i cr the series converges iff d > 2. 
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Now we want to study the scaling limit fj, — > fi cr + 0. Denote x = s\ + ... + sn- In the scaling limit 
one studies the exponents: mass (inverse correlation length) exponent is, susceptibility exponent 7, 
anomalous dimension 77, Hausdorf dimension dn- They are defined via the leading term behavior for 
small \i — n cr 

m(n) = lim lnG ( x ) ~ _ Mcr )^ 5 V G{x) « (p - ^ cr )~ 7 , 

|x|-t-oo \X\ ^ 

G{x) « |x|" <i+2 " ?J , NC(N; x) exp(-//iV) « 
We give now more detailed explanations. We have immediately that 



X ( M ) = V G(x) =Ve(-^) w = i 

V 4^ l-exp(-^ + /x c 



does not depend on the dimension of s and holds also for more general spins and interactions (however 
not always). Thus 7 = 1. 

If (p, — n cr )x <C 1 then G(x) Grw(x), the Green functions for the simple random walk. That is 
why 77 = 0. 

Let p( N \0,x) be the transition probabilities for the simple random walk on Z d and f(p,z) = 
^2n^2x p( N \0,x)e ipx z N ,z G C, - their generating function. Then 

/1 \ JV 1 

N d 
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which is the classical propagator for quantum relativistic free particle with mass \x — n cr . 
It can be also proved that 

d H = v~ x ,^ = v(2 - 77) 
In the second example different exponents can be obtained for A — > 00, fx — ► fi cr . 

2.2 Gravity as a queueing model 

One can construct a reversible dynamics with respect to which the measure Q(a) introduced above is 
invariant. It is a continuous time Markov chain. The state 

a = s\...sn 

is interpreted as a queue, where N is the length of the queue, S is the set of customer types and 
si + ... + sn is a generalized length of the queue (taking into account customer types and signs of 
jobs). This is a LIFO type queue (last in first out) and transitions consist of appending and deleting 
links on the right hand side (like arriving and service of customers in queueing theory). More exactly, 
for all a = s±...sn, with rate u 

Si...S N -> S1...SJV-I 

and with rate A 

si...s N -> Si...S N S N+1 

where all values of are equiprobable. Transitions from the empty queue are — > s± with rate 
A and equiprobable s±. 
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Lemma 1 If u > A then the process is ergodic and the distribution ^) is invariant with respect to 
this dynamics. 

To prove this, note that the restriction Nt = N(a(t)) of the process a(t) is also a Markov chain - 
a birth-death process, its stationary probabilities are 

Z~ l exp(- fj,N - N In 2), = l n ^ 

A 

For models 2 and 3 similar dynamics leaves the distributions invariant. Note that a system of two 
queues would correspond to two interacting particles etc. 



Supercritical case Now we see that they cases v < A have no sense in the equilibrium approach 
but in the dynamical picture they are no worse than the case v > A. One cannot write down equilibrium 
distribution for v < A, but at any time moment there exists some distribution Q(t) and its limiting 
properties as t — > oo could have interesting properties. Denote o~(t) the string at time t. 

Theorem 1 For A > v we have N(a(t)) — > oo a.s. Moreover there exist limiting local correlation 
functions (not too close to the ends of the string) which define a translation invariant Gibbs field. For 
example 

lim lim P(sk(t) = i) — > pi 

In fact for the first example this Gibbs field is a Bernoulli sequence in all three regions: v > A, v = 
A, v < A, see [45]. 



Critical case and scaling limit There are two possibilities to consider the critical case. 

The first one is to consider the dynamics for critical parameter values. The properties of this 
dynamics define the critical exponents. There are results for sufficiently general transitions: given two 
positive functions v(s), A(s, s'), define the transition rates as 

v(si...SN -> Sl...SjV-l) = Ksjv),A(si...Sjv -> Sl...SjV+l) = A(sjV,SAT+l) 

thus depending on the right symbols. Assume that the functions u(s),X(s,s') are such that the 
Markov chain is null-recurrent, see the conditions in [47]. Let S be finite with values a\, Let 
n(t) = (ni(t), ...,nk(t)),ni(t) is the number of symbols etj in the string a(t). Then 

Theorem 2 The central limit theorem holds for the random vector n(t), that is the following limit as 
t — * oo exists in distribution 

n(t) 

— = > \w\ c 

Vi 

where w has the standard gaussian distribution and c is a constant vector. 



This gives the same canonical exponents. Note that for the reversible case the proof reduces to the 
reflected random walks. The proof for non-reversible dynamics is more involved: for finite S see it in 
[47|. For compact S it should be similar. For non-compact it would be interesting to find examples 
with non-Gaussian limiting distribution. 

The second approach corresponds to the scaling limit in equilibrium case. In dynamics the param- 
eters are scaled together with time t, the parameters tend to the critical line v = A and v — A is scaled 
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as t 2 . In such dynamics the scaling limit corresponds to the diffusion approximation in queueing 
theory. One gets the Brownian motion with drift for the dynamics of x under the following scaling 

t = tN, x = rVN, v - A = iV~i 

The drift defines the mass gap in the spectrum of the infinitesimal generator of the corresponding 
diffusion process. The proofs here can be obtained by the application of the techniques known for the 
critical case. 

Random grammars We considered the dynamics, that is called right linear grammar (not 
necessarily context free) in the computer science terminology. Now we shall speak about more general 
dynamics when transitions can occur at any place of the string, not only in its right end. 

For the first example one can construct the following reversible Markov chain, leaving invariant 
the distribution, that appears to be a context-free random grammar (see |^3|]). Each symbol of the 
string is deleted with rate v and for each i = 0, 1, ...,n we insert a new symbol between symbols Sj 
and Sj+i (where for i = we put it before s\, and for i = n - after s n ) of the string s\...s n with rate 
A. Appended symbol with probability ^ will have one of 2d coordinate vectors e. To prove it note 
that this this dynamics restricted to Z + , the set of path lengths, is also Markov. It is in fact a birth 
and death process on Z + with jump rates qi.i+i = A(i — l),<ft$_i = vi,qoi = A. Then its stationary 
probabilities are (if A < v) 

A Ai...A n _i v 

7Tfc = ~ L exp(— fik), /i = m — 

u\...v n A 

For two other examples the dynamics (not context free) can also be constructed, we shall do it in 
another paper in more general cases. 



3 Spin Complexes 
3.1 Cell structures 

Here we present the minimum of basic definitions concerning cell structures. 

A complex is obtained by gluing together its elementary constituents - cells, like the matter consists 
of molecules. One should be very careful in defining the rules of gluing and the arising probability 
distributions. On the other hand it seems doubtful that some type of cellular structure has some a 
priori advantages in front of others. There are no definite physical reasons to prefer one cell structure 
or gluing rule etc., over another. Thus various possibilities should be studied to see what universal 
laws they share. In this paper paper we shall encounter two universal classes, one of them is popular 
in physics now. Moreover, having some flexibility in choosing a cell structure one can gain more 
simplicity in the probabilistic description and even get solvable models. 

3.1.1 Abstract complexes 

A (labelled) complex T is a set of elements called cells, there is a function dim A on T, the dimension 
of the cell A, taking values 0, 1, 2, .... The dimension of T is dimT = sup^ dim A. Let C T be the 
set of cells of dimension d. For each cell A G T^, d > 0, is defined a subset dA C U^Tq Tj, the boundary 
of A. Subcomplex T' of T is a subset of T such that if A S V then dA C T' . 

Isomorphism of two complexes is one-to-one mapping respecting dimension and boundaries. Equiv- 
alence classes of complexes with respect to these isomorphisms are called unlabelled complexes. 
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The star St(A) of the cell A is the subcomplex containing A and all cells B such that either A E dB 
or B E dA or dBndA^ 0. 

Note that complexes T can be considered as particular cases of spin graphs (G = G(C), s), seep4|]. 
The correspondence can be constructed in different ways. For example, let the vertices i of G corre- 
spond to cells of C, the function s(i) is the dimension of the corresponding simplex. Links are defined 
by the incidence matrix: two vertices A and A' of G are connected by a link iff A' E dA. 

Labelled spin complex is a pair (G,s) where G is a complex and s : C(G) — > S is a function on 
the set C(G) of cells of G with values in some spin space S. Isomorphism of two spin complexes 
is an isomorphism of the complexes respecting spins. The equivalence classes are called (unlabelled) 
spin complexes. Unless otherwise stated we consider only functions s defined on the cells of maximal 
dimension; by dualisation it is often equivalent to functions s restricted to vertices. 

3.1.2 Topological complexes 

There many topological incarnations of abstract complexes. In each of them a cell is represented by 
an open disk. A CW-complex is a topological space K which is defined by the inductive construction 
of its (i-dimensional skeletons K^. Let Kq = Tq be a disconnected set of points (vertices) - cells of 
dimension 0. In general, is obtained from Kd—i as follows. Each cell A of dimension d is identified 
with an open (i-dimensional disk Da and some continuous (attaching) map 4>a '• 9Da — ► K^-i is fixed. 
Then is the factor space of the union of K^-i and L>A:dimA=dDA via identifications of x E Da with 
(/>a(x). 

For example, K\ is a graph with vertices, zero-dimensional cells, and links (edges), one-dimensional 
cells. Link A is a loop if the boundary of A is mapped to one vertex. Often some restrictions on the 
attaching maps are imposed. Here we restrict ourselves to the case dimT = 2 and for all A, dim A = 2, 
the boundary 8Da is the union of some cells B, dimi? = 0, 1 (in some books, see for example Q, CW- 
complexes are defined as already satifying this restriction) . With such CW-complex one can associate 
an abstract complex with dA = {B E Kq U K\ : B C 4>a(Da)}- 

We get the class of simplicial complexes (where the cells are called simplices) if for each 2-cell A its 
boundary has 3 one-dimensional cells and the set dA n Kq of vertices uniquely defines A. Any graph 
without multiple edges, no loops is a simplicial complex. 

3.1.3 Cell surfaces 

In the paper we consider different classes of (two-dimensional) complexes. The class can be defined 
either by imposing further restrictions on the class of complexes defined above or by some constructive 
procedures to get all complexes in this class. Anyway such classes are a particular case of a language 
defined by some substitutions in a graph grammar, see |^3|, [P]] . 

The following restrictions hold for all complexes in this paper: complex is a (closed compact) 
surface. Pseudosurface (closed compact) is a topological space isomorphic to a finite 2-dimensional 
simplicial complex with the following property: each link is contained in the boundary of exactly two 
faces (two-dimensional cells). A surface has an additional property that the neighbourhood of each 
vertex is homeomorphic to a disk. 

This is the list of all compact closed (without holes) 2-dimensional surfaces. Orientable surfaces 
are just S p ,p = 0, 1,2, - sphere with p handles. Nonorientable surfaces are Pi (projective plane), 
Pi (Klein bottle), Pfc,..., - sphere in which k holes are cut and to each of them a Moebius band 
(crosscup) is attached along its boundary. 

In this case K\ is a graph homeomorphically imbedded to the surface S. Such complexes are 
studied in the topological graph theory (see ||) and in combinatorics, where topological complexes 
are called maps. Surface with holes is obtained from a closed surface by cutting out finite number of 
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disks with non-intersecting boundaries. If the surface has a boundary then the boundary belongs to 
K x . 

Isomorphism of maps is an isomorphism of abstract complexes. In other words, two maps are 
called isomorphic if there is a homeomorphism of S such that vertices map onto vertices, edges on 
edges, cells on cells. 

A map B is a subdivision of the map A if the graph K\{A) of A is a subgraph of the graph K\(B) 
of B. By Hauptvermutung if two topological complexes are homeomorphic as topological spaces there 
exist their subdivisions isomorphic as abstract complexes. 

If the surface is closed the Euler characteristics of the complex is defined as 

X=V-L+F 

where F is the number of faces, V- number of vertices, L - number of links. It does not depend on 
the complex but only on the surface itself: for orientable surfaces the Euler characteristics x = 2 — 2p 
where p is the genus (number of handles), for nonorientable surfaces x = 2 — n where n is the number 
of crosscups. 

We shall use in fact only the following 4 classes. 

Arbitrary maps This is the class we have just defined. No further restrictions are imposed. Sim- 
plest examples are a vertex inside the sphere (vertex map), an edge with two vertices inside the sphere 
- edge map. 

Smooth Cell Surfaces Smooth cell surface (see Q) is a compact connected smooth two- 
dimensional manifold M with finite number of closed subsets (cells) Fi such that: 

1. UFi = M; 

2. for each % there exists a one-to-one smooth mapping /, of F{ onto a polygon with > 3 faces; 

3. for i ^ j either Fj n Fj = or fi(Fi n Fj) is an edge or a vertex of the corresponding polygon. 

Triangulations This is a smooth cell surface with all rij = 3. The set V of vertices is called a cut if 
there are two subgraphs Gi, G2 such that G\ U G2 = G, G\ f] G2 = V. Disk-triangulation is a smooth 
cell surface, homeomorphic to the sphere, where there is one distinguished (that will be the outer face) 
face and for all other faces nj = 3, and there is no cuts with one vertex. Then it can be considered as 
the triangulation of the disk (sphere with a hole) . For triangulations the absence of cuts is equivalent 
to the absence of loops. 

Simplicial complexes These are triangulations without multiple edges, where moreover every three 
edges define not more than one cell. Note that a triangle (cycle of length 3) having inside and outside 
at least one vertex, is not considered cell. 

Convex polyhedra Quantizing smooth via piecewise linear structures is possible because the convex 
polyhedra have combinatorial counterparts. For example, convex polyhedra can be considered as 
maps with p = 0. There is a pure combinatorial characterization of maps corresponding to convex 
polyhedra. If a triangulation has no loops and no multiple edges, then, if L > 4, it corresponds (by 
Steinitz-Rademacher theorem), to a convex polyhedron. 
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Figure 2: Alexander moves 



3.1.4 Local observer (root) 

Labels in complexes are not necessarily given explicitly but the complex is considered to be labelled 
if the set V is claimed to be fixed. Labels are useful for fixing coordinate system in the space but are 
superfluous for the geometry and topology. There is a very convenient way to avoid the superfluous 
labelling but at the same time giving some algorithmic way to get a complete coordinatization. Root 
(local observer) in a (labelled) complex is an array (f,l,v) where / is a two-dimensional cell, I = 
its edge, v - vertex of /. Isomorphism of two complexes with roots is an isomorphism of complexes 
respecting the roots. Rooted map (rooted complex, complex with a local observer) of class A is an 
equivalence class of isomorphisms of complexes with a root in the class A of complexes. Assume that 
the rooted edge is directed from the rooted vertex. For disk triangulations we agree that one (the 
outer) face i is rooted, it is possible that m / 3. 

Lemma 2 The automorphism group of any rooted map is trivial. 

This is easily proved by induction on the number of cells by subsequent extending the automor- 
phism from the rooted face to its neighbors. 

3.1.5 Moves 

Graph grammars corresponding to transformations (substitutions here are called moves) of complexes 
were studied very little. In the next section we shall consider Tutte moves, see Fig. ^, which consist in 
appending an edge between two vertices of a cell or joining together two disjoint graphs by identifying 
two of their vertices. In topology subdivisions played always a big role. There are two papers (see 
0, ||]) where some moves are studied in detail. 

Let A be the commutative associative algebra over Z<i (simplicial chains over Z-i) generated by the 
symbols of some (countable) alphabet L with commutation relations sf = 0,SiSj = s,Sj. Thus it is a 
linear span generated by the strings (simplices) a = si...s n . Define the boundary operator as a linear 
operator d : A — > A such that 

8a = 

where the sum runs over all subsets of a with the number of elements |a| — 1. We shall consider here 
only two dimensional complexes. 

There are other linear operators in this algebra (called Alexander moves) Sij ]X ,i, j, x £ L. They 
are defined as follows. Let a = ij/3, then 

Si,j-, x ij(3 = x(i + j)(3 

Next example: Gross- Varsted moves. It is proved in (|| that each Alexandre move can be obtained 
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GV-move 1 




Figure 3: Gross- Varsted moves 



via Gross- Varsted moves and vice versa. We say that a set of moves is irreducible in the class A of 
complexes if for each pair T\ , Ti of complexes from A there is a sequence of moves giving T2 from 
T\ (in the physical literature the term ergodic is used in this cased, but we want to use the standard 
probabilistic terminology) . 

Theorem 3 In the class A of simplicial complexes the set of Alexander moves and the set of Gross- 
Varsted moves as well are irreducible. 

Proof see in f?], §. 
3.1.6 Automorphism group 

Let A be any of the five classes of complexes introduced above. 

Theorem 4 For most complexes with N two-dimensional cells from A the automorphism group is 
trivial, that is Anon £*^i (N) _^ g _> ^ where A(iV) (A non t r i V i a i(N)) is the set of all complexes 

with N two-dimensional cells from A (the same but with nontrivial automorphism group). 

Earlier Tutte remarked that it is very intuitive that almost all triangulations have no nontrivial 



automorphism. Many rigorous results appeared afterwards, see [18, 21, 22]. Proof for the case of 



disk-triangulations see in [23|. 



3.2 Metrics and Curvature 

The metric structure is defined once it is defined for each closed cell so that on the edges the lengths 
are compatible. There are two basic approaches for defining the metric structure: Dynamical Trian- 
gulations - when all edges have length one and Quantum Regge Calculus - when they are random. We 
shall use the first one. Then all cells with the equal number of edges are identical and on faces the 
metrics is standard. 

One can do it differently. Let first the graph K\ be embedded in the plane, the edges being smooth 
arcs. Define the metric structure on the graph K\ so that the edge lengths are all equal to a constant. 
Inside a cell with n edges we define the metric structure via some smooth one-to-one mapping of an 
equilateral polygon Q n with n edges onto this cell, so that the smoothness hold also in vicinity of each 
point on the edge. Then inside cells the curvature is zero. On edges also: this is shown on the figure 
in piecewise linear case. 

We shall define curvature R v at vertex v. As always the curvature is measured by parallel transport 
(Levi-Civita connection) of a vector (lying in the plane in piecewise linear situation) along a closed 
path: along the internal part of a triangle as on the Euclidean plane, through an edge - by unfolding 
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Zero curvature 



Positive curvature 



Figure 4: Curvature 



the two half planes separated by this edge to a plane. One sees immediately that only paths around 
vertices may give nonzero difference. Around the vertex v the angle between the initial and the 
transported vector is e v = 2tt — J2f Vfvi where ipj v is the angle of the simplex / at vertex v. Note 
that 

2*V - £v = E E Vfv = E E Vfv = * F 

V V f f V 

Using the Euler formula x = V ~ L + F one can get from this the Gauss-Bonnet formula 



^2 £ v = 2tt X 



for triangulations where L — "' 



2 



Gauss-Bonnet formula for smooth surfaces is J kdA = ^ J k(x)^/gdx = 2tt\ ■ Its relationship 
with the discrete case for a partition with nj-gons with areas A{ of the unit sphere (the area A of the 
triangle is A = ol\ + «2 + «3 — tt ) is given by the formula 



E ^ iAi ~ E 



^ Oij - {rii - 2)vr 



^2 ^ - n E ni + 27tF = 27rV ~ 2itE + 27rF 



Classical examples are: positive curvature - elliptic geometry (sphere, projective plane); zero 
curvature, Euclidean geometry (plane, torus, Klein bottle); negative curvature - hyperbolic geometry 
(all others). 

Now we shall show that the curvature Ri at vertex i is defined by the number qi of edges incident 
with i. Einstein-Hilbert action on the smooth manifold is 



/ 



(ciR + C2)y/gdx 

where R is the Gaussian curvature, g - metrics. It is known that j Rda = j R^fgdx = 47r%. Thus the 
discrete action should be (up to a constant) \p + fiN, where p is the genus and N is the number of 
triangles. We want to write down a discrete analog of this action with a discrete curvature summing 
over vertices instead of summing over triangles. Assume all triangles to be equilateral and scale 
their area to 1. Thus each vertex gets area | from each incident triangle, thus ^ in total. Then 
f Rda ~ J2i an d the formula 

Ef/"' 

i 

holds if only we put the curvature at the vertex i equal Ri = 2tt^-^. 
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4 Equilibrium planar pure gravity 



There are two kind of techniques used in the two-dimensional gravity. Historically the first one is 
the combinatorial approach, that was initialized by Tutte and continued (without any mention of 
physics) by many researchers, the papers are published in journals on combinatorics. The second one 
is Random Matrix Theory (RMT) approach, that was originated in physics itself. Calculations in the 
second approach are very persuasive but the arguments are not completely rigorous. As far as I know, 
no explicit connections between these approaches were established. We use the first approach and give 
a short review of the latter approach. 

4.1 Definitions and combinatorial approach 

Let A be some class of complexes (for example defined in the previous section) T, homeomorphic to 
the sphere, |T| = F(T) - number of cells of dimension 2 in T, C(N) = \ {T £ A : |T| = N}. The main 
example is the class of all triangulations of a sphere. 
The grand canonical ensemble is defined by 

P(T) = Z~ l exp(-/iF(T)), Z = £ exp(- / uF(T)) = ]T C(N) exp(-/xA) (2) 

TeA N 

In particular, the conditional distribution of T £ A with N fixed is uniform. Easy and general methods 
to estimate C(N) are useful sometimes, but can provide only bounds. 

Lemma 3 (exponential a priori bounds) c\~ff < C(N) < C27^, 1 < 7^ < 7 U < oo, Ci > 

Proof. Lower bound: this is quite trivial and can be proved in many ways. For example, take two 
following complexes homeomorphic to the ring with the same number of boundary edges from both 
sides. First one - alternating up and down triangles (that is standing on an edge and on the vertex 
correspondingly), second - two triangles up and two triangles down etc. These two kind of triangles 
can be glued sequentially one-after-one in all possible 2 n ways. 

The following method of proof of upper bounds works even in some more general situations. One 
can give an algorithm to construct all possible complexes with N cells of dimension two. Start with 
one cell. We enumerate its edges as 1,2,3. On each step we add not more than one cell to the 
boundary and enumerate new edges immediately after already used numbers. Now we describe the 
inductive construction. We take the edge with number one and make one of the 4 decisions: 1) not 
to add anymore triangles to this edge, 2) add to it a triangle having exactly two new edges, 3) add 
triangle to this edge and to the next edge on the boundary (in clockwise direction), 4) the same for 
counterclockwise direction. For each of 4 fc decision sequences uj let f(k,u>) be the number of edges 
after k steps, f(k,uj) < 3 + 2k. Moreover, if there are n triangles there cannot be more than 3n type 
1 decisions. 

One needs however exact asymptotics. All known examples exhibit the following asymptotic 
behavior 

C(N) ~ Cl N a c N (3) 

From (3) it follows 

Corollary 1 There exists < /x cr < oo such that for \x > n cr the series (||J converges. It diverges if 
fi < fi cr . If \i = fJ- cr then Z < oo iff a < — 1. 
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Thus, for the parameters fx < fi cr the distribution does not exist. However, the dynamics intro- 
duced later allows to consider such fx and for them local correlation functions make sense. 

No general results are known however. None of the constants c\ > 0, c > 1, is universal, but for 
all known examples a is. Universality of a is not at all simple intuitive fact. For example, predictions 
based on physical non-rigorous arguments (see, for example, [^]) failed to predict famous a = — \ in 
the planar case. 

Theorem 5 The asymptotics ^) holds for all four classes, defined in the previous section. Moreover 
a = — | in all cases. 

Proof. We shall prove it only for triangulations; other cases see in references cited in "Enumeration 
of two-dimensional maps" . In the similar way we shall define the distribution on the class Ao of rooted 
complexes 

P (T) = Z - 1 exp(- / uF(T)),Z = ^ exp(-/xF(T)) = £ C (N) exp(-fxN) 

TeA N 

where index zero means that we consider rooted complexes of class A. 
Lemma 4 For triangulations 

C(N) ~ (SN^CoiN) 

It follows from triviality of automorphism groups for most complexes (see theorem |]). Then we 
can take as a root any of N cells of dimension 2, choose one of its edges and orient it in 2 ways. 

Denote C(N,m) the number of disk-triangulations where the outer face has m edges, Co(N,m) - 
where the outer face is moreover rooted. The following result is similar but can be proved easier. 

Lemma 5 C(N,m) ~ m^ 1 Co(N,m) for large N and fixed m. 

Proof. Enumerate the edges of the boundary in a cyclic order: 1,2, ...,m. An automorphism <\> 
is uniquely defined, if j = (f)(1) is given. We shall show that almost all complexes do not have an 
automorphism such that j = <^>(1). 

To prove this we shall show that for each complex A having a nontrivial automorphism 4> we can 
subdivide the complex on two parts Ai U A<± = A where each cell belongs to only one part, such that 
(f>A\ = A2. This can be done by induction as follows. Take some boundary edge, take a triangle T 
with this edge and refer it to A\, then put (pT £ Ai- Each step of induction consists of taking one 
more triangle having common edge with already constructed part of A\ . Now we can modify Ai inside 
in a number ii(iV) of ways, bounded from below by some function u(N) — » 00 as N — > 00, uniformly in 
A. This can be done by choosing u(N) triangles in A2, not too close from each other, and modifying 
independently some neighborhood of each keeping the boundary of the neighborhood and the number 
of cells n in this neighborhood fixed. This is possible as C(n,b) > 1, where b is the number of edges 
on the boundary. Thus for given A\ the proportion of complexes with <pAi = A2 is small. 

We have proved that only "small" number of complexes have an automorphism (j) such that j = 
c/>(l). As m is fixed then multiplying this number on m gives again a "small" number. 

To prove the theorem we should prove that Cq(N) ~ C2N~^c N . The universal nature of (jH) is 
strongly supported by the fact that, for all such examples, the first positive singularity of the generating 
function Co(N)z N is an algebraic singularity, that gives the asymptotics (|3|) . 
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Remark 1 Assume an algebraic function y(x) is analytic at 0, has minimal positive singularity at 
point a > 0. We say that its leading exponent is b if there exist such bi and functions gi{x) analytic 
at x = a such that y(x) = Y^i=i 9i( x )( a ~ x ) bi + 9n+l( x )- Then we have the following expansion 



y( x ) = ^2 CnX ™' Cn ~ Cn ~ 



b - 1 a n 



In our case (for Cq(N)) b = §. One could also apply tauberian theorems in such situation. 

We give some examples where all constants in the asymptotics are known, see the same references. 

^3 
%■ 
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First example is the class of triangulations defined above. Here Cq(N) ~ 72-^ 2 c N ,c = 3-i/|. For 



convex polyhedra we have Cq(N) ~ 73 A" 2 c^ ,ci = For simplicial triangulations Cq(N) ~ 

ry 4 N~2c N ,c = ^p. Many other examples can be given; it is interesting however to understand the 
general underlying mechanism. 



Tutte |13| |l4j, |15[ has begun to study the asymptotics for C(N,m) and developed a beautiful and 



efficient "quadratic" method. Afterwards many authors contributed by developing the method itself 



and obtaining asymptotics for various classes A (see review [19[ and more recent papers p 



The main idea of Tutte are the following recurrent equations for C(N, m), N = 0, 1, ...; m = 2, 3, ... 
C(N,m) = C(N-l,m + l)+ ^ C(N 1 ,m 1 )C{N 2 ,m 2 ),m > 3, N > 1 

N 1 +N2=N-l,m 1 +m,2=m+l 

C(0,2) = 1,(7(0,771) = 0,m > 2 

These equations are easily derived as follows from the following picture where the orientation of the 
rooted edge is marked by arrow, rooted vertex is the first vertex of the arrow, rooted face is to 
the right of the arrow (containing the north pole of the sphere), see Figure Take any rooted map 
with (N — l,m + l) and do Tutte move 1, take any ordered pair of rooted maps (Ni,mi), (^2,7772) 
and perform Tutte move 2. Any rooted map (N, m) can be uniquely obtained in this way. (0, 2) 
corresponds to the so called edge map with one edge only which is counted twice. 
If we introduce the generating function 

00 00 

u(x, y ) = j2 J2 c ( N ' m ) xN y m ~ 2 

N=0 m=2 

the following functional equation 

U(x, y) = U(x, y)xy~ 1 + U 2 (x, y)xy + 1 — xy~ 1 U(x, 0) 

holds. We shall deduce from this equation that U(x,0) is algebraic and compute its first singularity, 
below in this paper, in a bit more general setting. 



Green functions Consider a class A of complexes. Let A(mi, m^) be a class of complexes, 
defined with the same restrictions as A, homeomorphic to the sphere with k holes with 777-j edges on 
the boundaries of these holes. We assume also that these boundaries do not intersect each other. The 
Green functions are defined as follows 

Z(mi, ...,m k ) = exp(-/xF(T)) = } ] C(N, mi, mfc) exp(- < uiV) (4) 

TeA(mi,...,m k ) N 
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Z corresponds to the case k = 0. Rooted Green functions are denned similarly 

Z (mi,...,m k ) = ^2 exp(-/iF(T)) = ^C (A r ,mi,...,m fc )exp(-/xiV) 
TEA (mi,...,m k ) N 

where the index everywhere means that we consider complexes with a distinguished edge on the first 
boundary with mi edges, the local observer in the terminology of p3| , |44|]. One would like to have an 
expression for the Green functions in terms of the basic probabilities (as for Markov chains). 

Green functions are associated with the derivatives = — x (m)> that is the factorial moments 
of .V. 

Lemma 6 The partition function and its two first derivatives are finite for \i = [i cr and for n > 2 we 
have as fj, — > fi cr + 

— — c(n) (M-Mcr) 



Proof. We shall see later that Z is an algebraic function of z = e~^ and has the principal singularity 
C(l — j^)~ a ~ 1 at the point z = e~ Mcr \ In the vicinity of z$ we have 1 — f - = 1 — e~^~ Mcr ) ~ /i — /i cr . 
This is in a good agreement with the following simple intuitive counting argument. 

Lemma 7 For fixed k,m\, ...,rrik there exist constants < c\ < C2 < oo such that 

dN^CoiN) < Co(N, mi ,...,m k ) < c 2 N k ~ 1 C (N) 

Proof. Take first k = 1 and prove the upper bound. Take some complex A £ Ao(m) with N faces 
and glue up the hole with some complex B € Aq(to) with r faces where r depends only on m. We 
shall get some complex C(A) £ Aq with N + r faces. For given A and C G Ao with JV + r faces we 
shall get not more than p complexes C(A) = C where p depends only on m. In fact, for any C the 
number of subcomplexes with r faces from Ao(m) having the same root is bounded by Cq(t, m). That 
is why Co(N,m) < cCq{N + r),r = r(m). The lower bound can be proved similarly. For k > 1 the 
proof is similar but one should first choose k — 1 faces along which paths with ui2 , • • • , rn^ edges will 
pass. This will give the factor N k ~ l . This can be done by induction in k. 

4.1.1 Uniform asymptotics 

Two questions arise: what is the asymptotics of C(N, m) if both N, m tend to infinity and what is the 
asymptotics of other global variables, such as the number of vertices etc. ? We shall see that these 
two questions are related. 

V,L,F are well-defined random variables in the grand canonical ensemble and one could would like 
to have their joint distribution. In general only two of them are independent due to the Euler formula 
V — L + F = 2. For triangulations, where each face has 3 incident edges, we have only one independent 
variable as L = y. For the class of all rooted maps, where two variables are independent, we have 
the following lemma. 

Lemma 8 Let E(V \ F = N) be the conditional mean number of vertices if the number of faces is N. 
Then 

E(V \ F = N) ~cN 

for some c > 0. 
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As it follows from the formula on p. 157 of [12| the number of rooted maps with N + 1 faces and 
m + 1 vertices is 

1 m 1 

c(N,m) = _ i C 2N+m-2 2 m -l C 2m+N-2 

Thus 

E(V | N) = ^mmc(N,m) 
E m c(N, m) 

is defined by the maximum in m = aN of ln(C^ +m _ 2 C2 m+ jv_2) by large deviation asymptotics. 
Consider now one-particle Green functions. 

Lemma 9 The following series 

J2e-^ N - um C{N,m) 

N,m 

converges above some nondecreasing function see Figure^ Thus the series Y^n m G~ flN C(N, m) 

diverges. 

Proof. It is quite obvious because the series has all coefficients positive. 



, M- 

Figure 5: Critical curve 

Thus we have a family of distributions P{fi,u). It is of interest to study the asymptotics and 
exponents when /i = av, where a, < a < oo is fixed. 

The explicit formula (see [jl2| ) for the number Co (N, m) of triangulations with a distinguished edge 
on the boundary (rooted triangulations) is 

_ 2J +2 (2m + 3j - l)!(2m - 3)! 
G ° (iV ' m) ~ (j- + l)!(2m + 2j)!((m-2)!)2 (5) 

where N = m + 2j is the number of inner cells. Cq(N, m) = if N — m is odd. 
Then as j, m — ► oo 



c » (iV -'")'- <:o "' i *j(ra) 2,+2 *" c -«, 



in particular for j ~ (3m , (3 > 



C (N,m) ~ C onst^ +2m -^^ 2+W)ml °^ 
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Thus for all < (3 < oo the exponent is a = —2. For fixed m the exponent does not depend on m 
and is a = — | and moreover 

C (N,m) ~ ^(mJJV-tc^ 
We have also <p(m) ~ Crn^c™ as m — > oo. 

4.2 RMT approach 

Random Matrix Model is the following probability distribution p on the set of selfadjoint nxn-matrices 
<fi = {4>ij) with the density 

^ = Z- 1 eM-H^) - tr{V)) 

where V = Yl a k<ft k is a polynomial of <p bounded from below, v is the Lebesgue measure on real 
n 2 -dimensional space of vectors {4> u , Re cj)^, Im^-, i < j). It can be written also as 

^ = Z^eM-HV)) 
"Mo 

where p is the Gaussian measure. It is easy to see that p has covariances {4>ij,4>*ki) = ( ( fiij,4'ik) = 
hSikSji. Note that for mere existence of the probability measure p one needs that the senior coefficient 
a p of V were positive and p were even. In this case there exists a well-developed probability theory of 
such models, which we shall not review here, see p7| . 

The fundamental connection (originated from t'Hooft) between RM model and two-dimensional 
complexes is provided by the formal series 

lo ^o = £ ^ < tr(V), ...,tr(V) >= E ^ E J (^) 

k=l ' ' D k 

where Dj, is the sum of all connected diagrams with k vertices. Take for example V = a^cfi 4 . Then 
each diagram has labelled vertices 1, ...,k, each vertex has labelled thick legs 1,2,3,4, corresponding 
to the product (frij&jk&kl&lv Each thick leg can be seen as a narrow strip with two sides, each side is 
marked with a matrix index. Dividing by 4 fc we eliminate the numbering of the four legs leaving them 
however cyclically ordered. After coupling legs and their sides (note that coupled sides have the same 
index and, as each vertex have two sides with the same index, we get index loops) and summing over 
indices we get a factor n N where N is the number of index loops. After this we are left with 

^ (~ Q 4) fc ^ h L{D ^n N{Dk) 

For each graph D choose the minimal cell embedding f(D) of D in a compact orientable surface 
of genus p (topological graph theory Q). Assume clockwise order of legs. It has k vertices, 2k edges 
and N faces. Putting a = —04 and using Euler formula k = N + 2p — 2 we have 



D k N,p 



= E C(N, p)(ah 2 n) N (ah 2 ) 2p - 2 = a 2 exp(-pN - up) 

N,p 

with p = — ln(ah 2 n), v = —2 ln(ah 2 ). 
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The calculations in RMM can be done only for n — > oo, thus to get finite \x one should scale as 
ah 2 = ^, b = e~^. In the limit n — > oo we have — > oo and only p = survives giving thus only plane 
imbeddings. The limit lim M ^ cr lim n _+ 00 is called the simple scaling limit. It was proved (see [Q) 
that a = — | in this case showing again the stability of this exponent. 

There are important points in this approach which should be mentioned: 

• In case V = a^ip the order of all vertices equals 4, this is some restriction on the class of maps; 

• Automorphism group of our labelled diagram factores in two factors. The first one C(D) related 
to the permutation of vertices, and second one Ci(D) related to the permutation of legs in each 
vertex. Almost all diagrams have the first factor trivial, but for some of them C V {D) > 1. We 
can sum over nonlabelled diagrams then but each unlabelled diagram will have a factor 

V Q(D) 
^ C V D) 

This means that the counting does not coincide with the natural counting used in the combina- 
torial approach; 

• We should fix also / somehow: normally one chooses embedding f(D) to the minimal possible p. 
But anyway not all possible triangulations are taken into account because a given graph can be 
embedded to surfaces with different p. This gives one more reason that the counting rule does 
NOT coincide with natural counting where all maps from some fixed class are counted exactly 
once. But this should not be taken seriously: anyway this counting is no worse and no better 
than others. 

• There appears a contradiction if one wants to get probability distributions simultaneously for the 
matrix model itself and for graph embeddings. We have probability distribution for the matrix 
model if a > 0, but the probability distribution on the diagrams is achieved only if a < 0. Thus 
one should always perform analytic continuation from a > to a < 0. The free energy for the 
scaling mentioned above can be rigorously calculated but the complete argument leading to the 
graph counting is still lacking. 

• There are other pure gravity models treated with this approach: more general pure gravity 
model counts the number n(q,T) of vertices v with q v = q: 

TeA<j>2 

where t q are the parameters, see |^|. 

5 Linear boundary dynamics 

The probability distribution 

P(T) = Z- 1 exp(-pF(T)) 

on some set A of complexes is invariant with respect to the following simple Markov process. Let at 
time t the triangulation be T{t). The process is defined by the following infinitesimal transition rates. 
With rate A+(iV), N = \T\, at time t we destroy T, add one more cell and glue anew all cells randomly 
together, that is if N = F(T(t)) then we choose uniformly T(t + 0) among complexes of the class A 
with N + 1 cells. With rate A_(A r ) we do random choice of a complex with — 1 cells. 
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Figure 6: Dynamics on the boundary 



What dependence on N can be ? If A+(T) = bf(N), A_(T) = df(N— 1) for some positive function 
f(N) then the probability distribution P(T) is an invariant distribution with respect to this process. 
Proof consists of the remark that the induced process on N is a reversible Markov chain: a birth and 
death process on Z + with jump rates q^i+i = &/(*), = df(i — l),<zo,i = 1- 

The simplest way of Monte-Carlo simulation is to take sufficiently large N and simulate uniform 
distribution, but it is impossible to find the exponent in this way. One should compare different N 
and this can be done via such a process. Apart from this such dynamics is of no interest, it is not 
constructive, especially in higher dimensions. In the rest of this paper we shall study local dynamics. 
We start with a simplest local dynamics of two-dimensional planar complexes. The distribution 
appears not to be invariant with respect to the first model dynamics. Thus, there could be two 
possibilities: either it will nevertheless give the same exponents for the invariant measure or its 
invariant measure belongs to another universality class (being however irreducible and ergodic). We 
shall show that the second one holds. 

5.1 Local Pure Growth 

We consider smooth cell surfaces and assume the cells be triangles. One starts with one triangle and 
each step consists in attaching a new triangle on the boundary. There are two kinds of attachment 
(see Figure 0) : to one or to two edges with the same vertex: for any edge on the boundary we attach 
to it a triangle with rate Ai. For any pair of neighboring edges on the boundary we attach to them a 
triangle with rate A2. At any time the complex is homeomorphic to a closed two dimensional disk and 
its boundary - to a circle. We assume that the initial state is the only triangle and that if the number 
of edges m on the boundary is equal to 3 then only Ai-transitions are possible. We can consider the 
states with m = 3 as giving a triangulation of the sphere itself (all other states as disk-triangulations), 
the outside of the triangle being the cell containing the north pole on the sphere. One can interpret 
it as the closing up of the hole in the sphere (the external part of the complex) . 

Remark 2 It is important to note that one could consider two other variants of this dynamics. First 
one is when we consider equivalence classes of cell surfaces. Then transition rates would be Ajm',i = 
1,2, instead o/Ajm, where m' is m divided by the number of automorphisms of the disk triangulation. 
Second, we shall use its analog later in more complicated situations, is that there is a distinguished 
(rooted) edge on the boundary and transitions can occur only if they touch this edge. 

There are 3 cases with quite different behavior of this Markov process: sub-critical or ergodic, 
critical or null recurrent, supercritical or non-recurrent. For all these cases we shall study the behavior 
of local correlation functions and of the following global variables at time t: 

F(t) = F(T(t)),V(t) = V(T(t)),L(t) = L(T(t)),m(t) 

where L is the total number of edges and m is the number of edges on the boundary. 
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Subcritical case By definition it is the case when A2 > Ai. 
Theorem 6 If A2 > Ai then 



A 2 



lim P(m(t) = N) = vr(m = N) ~jv-«x, CN~ l exp(-vN), v = In _ 

Let T3 6e i/ie random number of jumps until first return to the state 3 and put P{T) = P(T(t%) = T) 
for any triangulation of the sphere with a distinguished face (outer face). Then 

/A - A x N 
P(F(T) = N) ~ C 1 1 /l 



.Ai + A 2 , 

Proof. Note first that the length m(t) of the boundary is itself a Markov process: the evolution of 
the boundary can be seen as the simplest (context free) random grammar with the alphabet consisting 
of one symbol 1 (representing one edge) and with the substitutions 

Ai : 1 11, A 2 : 11 -» 1 

This process is obviously reduced to the branching process with one particle type where Ai is the 
birth rate, A 2 is the death rate. Denote this process m(t) - it is a continuous time Markov process 
states of which are the points of the lattice interval [3, 00). It has jumps ±1 and the corresponding 
rates X±k and X 2 k from the point k £ [3, 00]. Its initial state is 3. The process starts anew with the 
point 3. The stationary measure for the process m{t) is 

it(m = k) = ^\ =3 ^ ~ Ck~ l exp(-uk), v = — In — 

Remark 3 Note that for dynamics with a local observer the same considerations show that 

Y\ k ~ l A 

7r(m = k) = LLi = 3 1 ~ Cexp(-^) 

rii=4 A 2 

The exponents differ, as normal, by 1. 

Consider the (discrete time) jump process z n = m(t n ) for m(t), where t\ < i 2 < ... are the 
moments of jumps. Let f(s) = x 1 +\ 2 Ai-s 2 ) be the generating function for the transitions of the 
jump process and f n (s) be its iterates, let n(tu) be the first time when m(t n ) = 3. It is known that, 



see g|] 



P(n(u) = N) = (f N (0) - /at-i(O)) ~ Cm N ,m= ^—%L 

M + A 2 

This gives the proof. 

Supercritical case If Ai > A 2 then the boundary has exponential growth. The array (V, L, F, m) 
behaves like a degenerate branching process with four particle types and the following rates 

Aim : V ^ V + 1, L ^ L + 2, F ^ F + l,m ^ m + I 

X 2 m :V ->V,L-> L + 2,F -> F + l,m-+m-l 
Then. for some random variable ^ 

m(t) ~ ^ exp(Ai - X 2 )t, F(t) ~ c\£ exp(Ai - A 2 )t, V(t) ~ c 2 ^ exp(Ai - X 2 )t, 

L(t) ~ c 3 ^exp(Ai - X 2 )t 

taking into account the Euler formula V — L + F = 1 + m and L = the proof being the same as 



for similar statements in [44] 



23 



Critical case and the exponents If A2 = Ai then the process m(t) is null-recurrent. One 
cannot speak about its stationary probabilities but there exists an infinite stationary measure. It has 
however the same exponent a = — 1. 

We could consider a different analog of the stationary probabilities: we consider P(T(r) = N), 
where T(t) is the number of cells at the final moment r (when the hole closes up). 

Lemma 10 Let Ai = A2 = A. Then 



P(T(r 3 ) = N) 




It is well-known, | |46| |. 

Remark 4 The exponent is slightly different from the equilibrium case: a = —2 7^ — f • But also the 
following phenomenon occurs. Exponents in ergodic and null-recurrent case are different: a = and 
a = —2 correspondingly. That is the partition function for critical case converges and the limit of the 
partition functions as \i — > fi cr is infinite. 

Local correlation functions Even if the exponents are different from the physical theory it would 
be interesting to study local correlation functions, they define fluctuations of the curvature. 

Denote by i the i-th appeared vertex. Put qi(t) = if i did not appear before time t and otherwise 
put qi(t) equal to the number of triangles (or edges) incident to i. Let r(i) be the first time when the 
vertex i appeared. 

Theorem 7 There exist XkiYlXk = 1> such that for X\ > A2 

lim lim P(qi(r(i) + s) = k) — > Xk 

s — >oo i — >oo 

For some constants a,b > and any two vertices i < j such that the initial distance p(i,j) (that is 
the distance at time t(j) ) between them equals d 

\P(qi(t) = k,qj(t) =1) -XkXil < bexp(-ad) 

for all i and t sufficiently large. If X\ > A2 then 

lim lim P(qi(r(i) + s\m(r(i)) = m) = k) — » Xk 

s— >oo m— »oo 

Also the same exponential decay property holds. Moreover, all local correlation functions Xk are ana ' 
lytic functions of X\, A2 for all values of Ai, A2 > 0. 

Proof. For the vertex i put £j( s ) = qi(r(i) + s). Note that £j(0) = 2 and the number % will increase 
until both adjacent links to this vertex will not enter one new triangle, denote this random time <Tj, 
it has exponential distribution with the parameter A2. Thus Xk ^ s equal to the probability that the 
Poisson process with rate 2Ai + 2A2 will have exactly k — 3 jumps, that is appending triangles to one 
or two vertices from the left (or from the right) of i. In fact, take vertex i and its two edges. Until 
this moment from both sides new triangles in i appear. The crucial argument is that this process is 
independent of m as far as during this time interval m > 3 and independent on all events which do 
not touch vertex %. 

The second assertion of the theorem is proven quite similarly if we remark that the processes q%{t) 
and <7j(t) become dependent only when the distance between i,j becomes less than 4. For this to 
occur there should be many A2-events in-between % and j, which has exponentially small probability. 
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Define the mean curvature 

From the exponential decay it follows the central limit theorem for the scaled curvature 
Theorem 8 For any sequence of sets I of vertices 

T,ieiRi-W\ 

converges to Gaussian distribution as \I\ — ► oo. 
5.2 Reversible boundary processes 

The reasons why for this dynamics we did not get the desirable invariant measure on the complexes 
are rather delicate. We shall give now an intuitive explanation. One could expect a simple invariant 
measure for a process which is reversible. To get a reversible process we should add the possibility 
of deletion of faces. For example, let Ai = A2 = A and assume that each boundary triangle can be 
deleted with rate \i unless the resulting state does not belong to the class A of complexes. 

Lemma 11 This Markov process is reversible. 

Proof. Let Qj be complexes. We shall consider closed paths ati, 0:2, ct n , a\ where each ccj+i is 
obtained from aij by appending or deleting a triangle on the boundary. The number of appending in 
such closed path cti, ct2, ct n , ct\ should be equal to the number of deletions, thus it is ^. Denote X a b 
the rate of transitions from complex a to complex j3. Thus 

^ai02 , "^«i i ai ^aia n -- -^012011 (-^A*) 2 

There are however complexes with arbitrary N where no triangle can be deleted, because otherwise 
we get states, where two graphs intersect only in one vertex, which are not allowed. The set of such 
states has a sufficiently complicated nature and it is difficult to use standard procedure to get stationary 
probabilities via balance equations. However, we could split these graphs into two parts belonging to 
our class, then we can get splitting the complex into connected components. It is exactly the latter 
operation which brings us to nonlinear Markov processes. 

6 Nonlinear boundary dynamics 

In the preceding section we considered a Markov dynamics, the states were the complexes themselves. 
This dynamics was a local dynamics: the changes could occur at any point of the boundary. We saw 
that this did not give us exponents accepted in the physical literature. Here we will construct dynamics 
giving exactly the exponent a = — |. This dynamics appears not a Markov process. One can think 
about infinite number of universes interacting with each other, but not like an infinite particle system. 
Simplest pairwise interaction (gluing) of two universes will give us quadratic functional equations and, 
as a result, the necessary exponents. 
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6.1 Quadratic quasi processes 

Let A be the set of probability measures on some space S. We shall consider a class of transformations 
M : A — > A, generalization of Markov chains. These transformations are not generated by random 
maps S — ► S, they are nonlinear on A. 

Let a (denumerable) set S be given and let a Markov chain be defined on the state space S with 
transition probabilities p a p, Y^p-p^aPafi = 1> 110111 a to /3. It defines a linear transformation L 

9 = {9/3} Lq = q a p a p} 

a 

on A. Let also a probability kernel P((a,7) — > 0) : S x 5 — > S be given, -P(( a >7) ~* P) = 1- It is 
symmetric P((a,7) —>/?) = P((7, a) — > /?) and P((a,7) — ► 7) = 0. It can be deterministic. It defines 
the quadratic transformation Q on A 

q = Ua} ^Qq = {^2 WiPdPi 7) ")} 

/3,7 

Then taking a convex combination we have the transformations (formally cq + C\L + c 2 Q) on the class 
of measures on S 

q(a, t + l)= ri J2 *)P/3a + ^2 J] ?(/9, t)9(7. t)P((j3, 7) -» a) + (1 - n - r 2 )c (a) (6) 

for some probability measure Co (a) and nonnegative numbers r\ , r 2 such that < r\ + r 2 < 1 .We see 
that the total mass q{a) = 1 is conserved. 

This can be interpreted as follows. Consider a denumerable number of particles on S, q a is the mean 
number of particles at point a. With probability 77 each particle (independently of the others) makes 
a jump according to the probabilities p a p. This gives the linear transformation. With probability 
r 2 particles form pairs so that mean number of pairs (a, (3) is the product of means, then each pair 
(a, 7), independently of each other, gives birth to one particle at (3 with probability P(a,7 — > 
Also with probability 1 — 77 — r 2 one has an immigration with mean cq{o) of particles to a. We want 
to emphasize that there is no stochastic process here in its standard sense but only the transformation 
of measures. This seems to be related also to field theory of strings (second quantization of strings) 
but this physical theory does not have a mathematical status. 

Some probabilistic theory There is no theory of such quadratic quasi-processes and we have to 
discuss it here. We can rewrite it in more general terms 

T = Pl Ti + P2 T 2 + (1 - Pl - P2 )T 3 

Let k\ be the contraction coefficient for T±, that is for any probability measures fa, fa we have 
11^1(^1 — ^2)11 ^ &i 1 1 Mi — fa II- Let fe 2 (x) be the contraction coefficient for the stochastic matrix 
Pyz(x) = P((x,y) — > z) and k 2 = sup x k 2 (x). 

Theorem 9 Assume p\k\ + 2p 2 k 2 < 1 . If the number of states is finite then there is exactly one 
fixed point of T and the convergence to it is exponentially fast. If the number of states is countable 
the same assertion holds under the condition that there exists a (Lyapounov) function f(x) such that 
^2 x f(x) < 00 and if (i(x) < f{x) then also (Tfj,)(x) < f{x). 
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Proof. Take two probability measures v and fi = v + e. Then we have the following contraction 
property for = ||e|| 

p(Tfi, TV) = | Pl Tie + p 2 ^ Kx)e(y)P((x, y) -► z) + p 2 ^ u(x)e(y)P((x, y) z) 

< (pih + 2p 2 k 2 ) \\e\\ 

Then the first assertion of the theorem follows. To prove the second assertion note that by compactness 
there is a fixed point v in A = {/j, : /j,(x) < /(x)}and for each /j, the sequence T n ^ converges to v. 

Continuous time quasi-processes are defined similarly. Instead of probabilities p we introduce 
rates A: with rate Ai we do the linear transformation, with rate A2 the quadratic transformation, and 
immigration arrives with rate Ao- The equations for the stationary measure are the following 

(A x + A 2 + Ao)vr(a) = Ai ^ tt((3, t) Ppa + A 2 n(P, t)ir(-y, t)P((P, 7) -> a) + A c (a) 

0,1 

and can be reduced to the previous case. The time evolution is governed by the following equation 

= Ai t) - q(a, t))pp a + 



+A 2 ^2(q((3, i)g( 7 , t) - q(a, t))P((p, 7) a) + A (c (a) - q(a, t)) 

0,1 



6.2 Generating functions 

Now we consider the dynamics with the set S of all plane disk-triangulations with the root on the 
boundary and the projection of this dynamics onto where the points of Zl. are denoted a = (N, m), 
N is the number of faces (not outer) and m is the number of boundary edges. That is the measure of 
the point (N, m) is the sum of measures of the corresponding complexes. This projection will appear 
to be also a quadratic quasi-process. Introduce the generating function 

00 

U(x,y)= q(N,m)x N y m 

N,m=0 

We assume homogeneity: for all a, (3,a±,(3 1 ,^ such that a, (3, a + 7, P + "f, ai, Pi, cci + (3 1 +7 all belong 
to the quarter plane 

Pa,a+j = a 7 , P{(a, P) -> a + /? + 7) = P((ai, -> ai + P 1 + 7) = 6 7 

together with the following bounded jumps assumptions: p a p = if \a — P\ > d for some fixed integer 
d, P((a,P) — > 7) = if \a + P — j\ > d, co(a) / only for finite number of a. Introduce the 
generating functions 

A(x, y) = Y^ a ( n, m )x n y m , B(x, y) = Yj b (n,m)x n y m , C(x, y) = Y, c o( n , m)x n y m 

n,m n,m n,m 

Now we get a functional equation for the generating function of the stationary measure 
U = riU A(x, y) + r 2 U 2 B(x, y) + (1 — 77 - r 2 )c (x, y) + b.c. 
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Linear move 



+ 





Quadratic move 



Figure 7: Linear and quadratic moves 



where boundary terms can appear because there is no homogeneity in the vicinity of axes of Z 2 ^. If 
d = 1 then the boundary terms are linear combinations of functions U(x, 0) and U(0, y) and finite 
number of ir(N, m). 

Consider now the case where the jumps correspond to the moves shown on the Figure [?] 
Here the jumps are 

P(N,m),(N+l,m-l) = 1 

P((N 1 ,m 1 ),(N 2 ,m 2 )) -» (#i + N 2 + l,mi + m 2 - 1) = l,c ((0,2)) = 1 

To calculate boundary terms it is convenient to introduce a complex consisting of one edge with two 
vertices, it corresponds to N = 0, m = 2. Then cq is the unit measure on this complex. Gluing of 
two such complexes gives a triangle with N = l,m = 3. From the line m = 2 linear jumps are not 
possible, this is the only line where the homogeneity is destroyed. Thus the equation is 

U = nlJxy' 1 + r 2 U 2 xy' 1 + (1 - n - r 2 )y 2 - xy^ri ^ q(N, 2)x N y 2 

N 

One can do the following scaling in the main equation 

x — * xr^ 1 , q(N, m) — * aq(N, m) 
If we put ar 2 = 1, (3 = ^ 1 ~ ri ~ r " 2 ^ r2 then it becomes 

U = Uxy- 1 + XJ 2 xy- 1 + (3y 2 - xy' 1 ^ q(N, 2)x N y 2 

N 

We call this equation canonical. Note that the case /3 = 1 or r± = r2 j-^_~^ corresponds to the counting 
problem and was completely solved by Tutte. We only reproduce his analysis in a more general setting. 

If there exists a solution U with nonnegative coefficients q(N, m) of the canonical equation then 
the invariant measure q(N,m) for the dynamics will be q(N,m) = a~ 1 r^ q(N,m). 

We shall be interested only in the class E of invariant measures satisfying the following assumptions: 

1. (exponential bounds) q(N,m) < C N+m for some C > 0; 
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2. q(N, m) are nonnegative 



The main result is the following theorem. 

Theorem 10 A unique positive invariant measure exists for all r\,r%. It is finite iff || 2 ( 1 - r ^-^) r ^ < 
1 . The measure of the sets {(N, m)} for fixed m is finite iff 2 ( 1 < \ 

In fact we have to prove the following 

Lemma 12 For any (3 there exists a unique solution of the canonical equation with positive coefficients 
in the considered class of measures. The series q(N, m) < oo converges for each m iff j3 < ^ . 

The series Y^n m?(^' m ) converges iff%\J^p > 1- 

Proof of the Theorem. One can easily prove that for any initial measure the difference of measures 
for any two complexes with the same (N, m) tends to zero. Thus, any invariant measure has the 
property that all complexes with the same (N, m) have the same invariant measure. Then we can use 
the above scaling argument. The theorem follows. 

6.3 Analysis of the functional equation 

Here we prove the Lemma. Putting U = y 2 W and defining S m (x) by 



W(x,y) = J2 S m-2(x)y % 



m=2 

we can rewrite the functional equation 



W = p + xyW 2 + xy~ l {W - S), S(x) = S (x) = W(x, 0) (7) 

In this case there is only one boundary term S (this simplifies strongly). This is the generating 
function for the number of disk-triangulations with rooted face, having exactly 2 boundary edges. In 
the general case the equations would be more complicated. It demands joining together the quadratic 
method of Tutte and the methods developed by the author for random walks in a quarter plane. 
If we know Sq(x) as a formal series then all S m are defined recursively by 

yS = f3y + yxS 1 ,y 2 S 1 = y 2 xSl + y 2 xS 2l ... 

We partially follow the derivation in [12]. We rewrite the functional equation (|7|) in the form 

(2xU + x - y) 2 = Ax 2 y 2 S + (x - y) 2 - 4(3xy 3 = D (8) 

Consider the analytic set {(x, y) : 2xU + x — y = 0} in a small neighborhood of x = y = 0. Note that 
it is not empty, (0,0) belongs to this set and it defines a function y(x) = y(x) = x + 0(x 2 ) in the 
neighborhood of x = 0. We shall prove that y(x) and S(x) are algebraic functions. We have two 
equations valid at the points of this set 

dD 

D = 0, — = 

dy 

or 

4x 2 y 2 S(x) + (x- y) 2 - 4/3xy 3 = (9) 
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8x 2 yS(x) -(x-y)- I2[5xy 2 



from where we shall get both y(x) and S(x) 



x = y(l-2(3y 2 ),S 



13(1 ~ Wy 



(1 - 2/3y2)2 

The algebraic function y(x) satisfies the equation y 3 + py + q = with 

1 x 

Its discriminant 

A 



-4p 3 - 27g 2 



p 2 (--27x 2 ) 



is not a square in the field of rational functions. Then (see 
ramification points are x = db 



the Galois group is S3 and the 




2I0. Cardano (formal) solution of the cubic equation is 




Note that If I ^ 
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and thus there is no singularities inside the circle of radius x\ = y 2fp- 
Two terms with opposite signs give the cancellation of the lowest order singularity (x — xi)z and we 

3 

have thus the leading singularity (x — x\)^. We need the branch where y(0) = 0, as y(x\) > then it 
is known that if the discriminant is zero then y(x\) = Then iterating the equation 

x 

y 



1 - 2/V 

we get that the expansion of y(x) at x = has all coefficients positive. S is an algebraic function, 
analytic for \x\ < y^j- I n l&ct, S could have a pole for |x| < y^f-g only if 1 — 2(3y 2 (x) = but it 
would imply x = which is impossible. To visualize the expansion of U denote now yo(x) = y(x) and 
substitute 

P(l-3P 2 y 2 (x)) 



S 



,x = y (x)(l-2Py 2 (x)) 



into 



[\-2l3y 2 (x)f 

We get 

D = Ay 2 y 2 (l - 3(3 2 y 2 ) + (y (x)(l - 2f3y 2 {x)) - y) 2 - 4(3y 3 y (x)(l - 2f3y 2 (x)) 

= {y-yo) 2 (a(x) + b(x)y) 
as y = yo(x) is a double root of the main equation, and 

a(x) = (1 - 2Py 2 (x)) 2 ,b(x) = -4/3y (x)(l - 2(5y 2 (x)) 
Choosing minus sign we have then 

y-x \f~D -x + y n (x) , „ \/a(x) + b(x)v — 1 



U(x,y) 



-x + y Q (x) v / a(x) + b(x)y 

- [y - yo{x)) — 



2.T 



2x 



2x 2x 
that gives a legitimate expansion. 

For given x the convergence radius of U as the function of y is defined by zeros of \J a(x) + b(x)y 

or Jl-^y. As — increases on the interval [0, x\] then the convergence radius for x = x\ is 
R = Aa *i_ f . We have ,R = Ux v 
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6.4 Correlation functions 



Using the tree representation introduced below we prove that for most vertices the conditional distri- 
butions of the random variables q v converge to the unique limit as N — > oo. One should know how to 
specify a vertex v if they are not labelled. Below we give some way to do it, we discussed this problem 
in p3|, 44] in a wider extent. The proof is combinatorial but without use of analytic methods, it is 



sufficiently involved and we present it not in a completely formal way. 



6.4.1 Tree representation 

The generation process of maps, given by recurrent application of Tutte moves, will be represented as 
a planar tree. Moreover, this will give a one-to-one correspondence between maps and some class of 
planar trees. A planar tree has a root vertex and grows upwards (it is shown on Figure || by arrows). 
Denote Ao the class of rooted disk-triangulations. We shall denote vertices of trees by v and vertices 
of maps by w. 

Denote To the class of all planar trees characterized as follows. There can be 3 types of vertices: 
0, 1, 2 according to how many edges go upwards from this vertex. Vertices of type are also called 
end vertices. Denote rii the number of vertices of type i. Among 2-vertices there are vertices which 
are incident to one 0-vertex, let their number be H2q and which have two incident 0-vertices, their 
number is denoted by uqq. The only further restriction on this class of trees is the following. For 
any vertex v denote T v the tree consisting of the vertex v and all vertices above v. Denote rii(v) etc. 
- the corresponding numbers for the tree T v . Note that uq(v) = n2(v) + 1. Then the class To is 
characterized by the following restriction: uq(v) — ni(v) — 1 > for all vertices v of type 1. 

Lemma 13 There is a one-to-one correspondence between Aq and Tq. 

Show first that each map generates a planar tree in To, that is there is a function / from maps 
to planar trees. We prove this by induction on the number of faces. The map itself is represented 
by the root vertex of the tree. From the root vertex we draw upwards one edge in case of move 1 of 
Tutte, and two edges, in case of move 2, corresponding to splitting the map on two maps. We can 
distinguish the latter maps corresponding to the right or left vertices of the tree accordingly to the 
orientation of the rooted edge of the map. In fact, the rooted edges of the two maps, on which the 
map is split, have the same orientation. Thus one of then precedes the other one. The map with the 
preceding rooted edge we consider as corresponding to the left vertex (with smaller numbers) and the 
map with subsequent rooted edge - to the right vertex. Each step reduces the number of faces by 1, 
thus we come to induction hypothesis. If the map is the edge map (consisting of one edge), then the 
corresponding vertex is the end vertex of the tree. 

Let us note that N = n% + ri2, L = uq + n\ + n2, V = no + 1, m = n$ — n\ + 1. Here N, L, m, V 
refer to the map, and numbers nj,noo etc. - to the tree. To prove the latter equality note that each 
1-vertex and each 2- vertex (except those which are incident to 0-vertices), when passed downwards, 
diminish m on one, each of 77,20 vertices gives one more edge to m and each of noo vertices gives three 
more edges to m. Thus m = n2o + 3noo — n\ — (712 — uqq — n2o) = no — n± + 1- As m > 2 for each 
map, we have the restriction. 

Now show that each planar tree of class To generates a map - we shall show that / is one-to-one. 
Take a tree and proceed by induction from upper to lower vertices. All end vertices we declare to 
be edge maps. Take V end edges, enumerate them from left to right as the end vertices of the tree 
v = 1,2, V. These edges will give V + 1 vertices in the complex. We mark all end vertices. Each 
induction step we take a vertex v such that there are only marked vertices above it. If v is of type 
i then the induction step consists of Tutte move %. We mark v after this step. Then we proceed by 
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Figure 8: Planar trees and maps 

induction. Inversely, the map constructed in this way generates the tree from which we started. All 
maps are legitimate because m > 2. Lemma is proved. 

Contribution of the tree G(T) is defined as the product r^V^V^ 2 where rii = nj(T) is the number 
of i-type vertices. This is equal of course to the probability of the corresponding rooted map. 

Remark 5 Planar trees are in one-to-one correspondence with the parenthesis systems, which can 
be put the product a\...a n in non associative non commutative algebra, see Figure || However, the 
restrictions posed on Tq make this one- dimensional grammar more involved. 



- ((()()))(()()) 



Figure 9: Planar trees and parenthesis 
6.4.2 Local curvature 

Take some end vertex (or end edge) v and put l(y) = max(u, V — v). 

Theorem 11 If l(v) — > oo (necessarily V — > oo and N — > oo) then for the equilibrium distribution P 
there exists the limit 

oo 

\im.P{q v = k) = p(k), ' S ^p{k) = 1 

k=2 

Proof. One can imagine V edge maps, corresponding to the end vertices, put horizontally along 
the line and directed from left to right. Take some v and the left (that is the rooted) vertex w = w(v) 
of the corresponding edge map. We fix orientation of the boundary to be counterclockwise. Define 
the history of each vertex w(v). This history can be described in terms of maps and in terms of the 
tree, it will be convenient to use both descriptions. By definition it consists of several parts, that we 
call history parts, of the unique path from v to the root of the tree in the tree T in the downward 
direction. In the sequence of maps the history ends either when the vertex w(v) disappears from the 
boundary of the map (this means that q w will not change anymore) or in the root of the tree, if w is 
on the boundary of the final map. The first part of the history lasts until the vertex will be covered 
by an edge. q w may change only on such history parts. From the tree point of view the first part of 
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the path starts with the end vertex v, goes downwards vertex by vertex, and lasts until the first right 
2-vertex.We say that the 2-vertex is a right vertex if our vertex w(v) is in the right map. Next part 
of the history starts when w(v) becomes to belong to the rooted edge of the map and ends exactly as 
in the first part. 

If the vertex a of the tree is above v = w(v) then w belongs to the corresponding map f~ 1 (a). 
Denote a(w) the infimum of a such that a < v and w is on the boundary of / _1 (a). 

Denote b(w) the infimum of b such that b < v and w is still a rooted vertex in the corresponding 
complex for each bi on the interval v > b\ > 62 > ••• > b in the tree. The length L(v(w)) of the latter 
path in the tree is exactly the difference q w (j) — q w (k) where j is the first complex of the path (edge 
map), k - the last one. We shall prove that Py(L(v(w)) = n) < C exp(— 771) for some 7, C > 0. Note 
that we can use for this any invariant measure as this conditional probability is the same for fixed 
N, m (in fact we are interested only in m = 3). 

If from the tree one deletes all 1-vertices then the resulting tree without 1-vertices will be called a 
bare tree, It defines an equivalence class of dressed trees, each of them can be obtained by appending 
some number of 1-vertices to the bare tree. Each nonnegative measure on trees induces a measure on 
equivalence classes - bare trees. We start with bare trees. 

Case of bare trees It is the case when there are no 1-vertices at all, that is m = no + 1 = N. In 
this case all vertices are on the boundary, and once the vertex was covered it does not participate in 
the process anymore. It can be covered only when a left join occurs and until it was covered there can 
be k right joins which give k extra edges to q w . Any right join joins to w a complex corresponding to 
some tree covering an interval to the right of v. 

For example, let us estimate the probability P(k;ji,...,jk) that there are exactly k joins to the 
intervals of lengths ji, ...,jk covering the interval I r = {v + 1, v + R} before the left join covering 
the interval 1\ = {v — L, ...,v — 1}. We shall use the following Markov property. Call a vertex v' = v(I) 
/-separating if it covers exactly the interval /. Then the number of trees with the separating vertex v' 
is equal to t(I)t(V — I + v ') where t(I) is the number of trees on the interval I, because t(V — I + v') 
is also the number of factor trees with respect to the set of trees on I. Thus the probability that v' is 
/-separating can be estimated as 

t(I)t(V -I + v') . 3 
t( yj -~min(6,c7 2) 

as t(n + 1) = -^yyCV; n ~ cn~^2 n are well-known Catalan numbers and b < 1. It follows that the 
probability P v that the vertex v will have a right join earlier than a left join. It is clear that P v < d < 1. 
Then using the Markov property we shall get by induction for large r and some a < 1 

P{q v >r)<a r 

Thus we got the exponential estimates. The existence of the thermodynamic limit follows from the 
fact that the influence of the boundaries takes place also with probabilities less than a l( - v \ 

Case n\ > 0. We fix a bare tree and consider one auxiliary problem (urn problem) concerning 
the distribution of 1-vertices on the bare tree. The estimates are uniform in bare trees (equivalence 
classes) . 

Consider first the probability that the vertex w will get large value of q w due to 1-vertices until it 
will covered at the first time. We shall do such estimates separately for each history part and consider 
in detail only the first history part. Let v be the vertex where the first parts ends. Let T v be the tree 
over this vertex (with the root v) and f~ l {v) be the corresponding complex. 
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First consider the case when T v has only one 00- vertex, then 0- vertices join sequentially to already 
existing trees. It means that different ways to put 1-vertices to the bare tree can be identified with 
all possible ordered arrays of nonnegative integers a±,...,a n such that for all k = 1,2, ...,n we have 
Eti a>i < k - 2. 

We can formulate the following abstract urn model. Let we have n urns and m balls in these 
urns, aj - the number of balls in the urn i. Let c(n,m) be the number of arrays a±,...,a n such that 
Yli=i a i ^ k> YliLi a>i = m < n. Then we have the following recurrence 

m 

c(n, m) = c ( n — l,m — k) 
k=0 

or 

c(n, m) = c(n, m — 1) + c(n — 1, m — 1) 

from where it is not difficult to get asymptotics for the number c(n,m;i) = c(n — l,m — i) of arrays 
among c(n, m) such that there are exactly i balls in the last urn. Then we have an explicit formula 
for the generating function 

oo oo oo 

f(x,y)= £ E^^y ^-l^ E^ 

n=-oora=l yV ' fc=l 

the coefficients coincide with ours for n > m. We want to prove that c ( ra >" 1 ) < 5 f or some 6 < 1, 

— r c(n,m—l) ' 

and to find a method (not using generating functions) which could work in the general situation. For 
this we rewrite c(n, m) in terms of the number of paths starting at the line m = 1 and ending at the 
point (n,m). Steps of the paths are either (0, 1) or (1, 1). We have 

n 

c(n,m) = kL(k;n,m) 

k=m+l 

where L(k;n,m) - the number of paths from the point (k, 1) to the point (n, m), as c(k, 1) = fc. As 
L(/c; n, m — 1) = L(/c — 1; n — 1, m — 1), then ^-1^-1) = 1 + O(^). The result follows. 

Consider now the general case. Instead of the urn problem on the interval [l,n] we have ah urn 
problem on an arbitrary planar tree T under the conditions 

^2 a v < V(T V ), ^2a v =m 

where a v is the number of balls in the urn (vertex) v of the tree, V(T V ) is the number of vertices of 
the tree T v . Let c(T v ,m) be the number of such arrays on the tree T v with m balls. If for example 
from the vertex v only two edges go upwards to the vertices v(l),v(2), then 

m 

c(T v ,m) = ^2c(T v ,m;i),c(T v ,m;i) = ^ c(T v{1) ,m 1 )c{T v ^),m 2 ) 

i=0 m\+m2=m—i 

Then the argument is similar to the previous one. Note that c(T v , 1) = V(T V ). We want to compare 
c(T v ,m; i) and c(T v ,m; i + 1), for this we iterate the latter recurrent equation for c(T v ,m; i + 1) to the 
very end, that is we get the sum of terms B l s +1 , in each of them all factors equal c(T v > , 1) for some v'. 
The iteration process for c(T v ,m;i + 1) there corresponds the similar process for c(T v ,m;i), that is 
why to each term E>1 +1 there corresponds the term B' l s in the expansion for c(T v , m; i). In that term one 
of the factors is c(T v /,2) instead of the factor c(T v /, 1) in the term Bl +1 . Thus as before ?^F V '"^ > 2. 
From this bounds uniform in bare trees follow. The influence of the boundary is exponentially small. 
Similarly one can estimate other correlation functions, for example, the decay of correlations. 
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Theorem 12 Let V — ► oo and take two vertices Vi,V2 with l(vi) — > oo. T/ten 

Kfefe) - < cexp(-a|wi - v 2 |) 

7 Internal dynamics 

We considered above only a growth of the boundary, that was quite natural: many modern technologies 
follow this principle. But also another dynamics is possible where all cells (even inside the building) 
can evolve. We shall consider here some questions related to such dynamics. 

Note that Gross- Varsted moves can be used not only for simplicial complexes but for other classes 
as well, as it is seen from the picture. Consider GV-moves 1 and 2 and the inverse one to 2, consider 
the Markov chain with rates Xi,i = 1,2,// for these moves correspondingly. 

Thermodynamic limit of local processes If A 2 = /t = then V,L,N are invariants. Let 
A' C A(N, L) be an irreducible component of the set of (nonequivalent) complexes with given TV" and 
L and C(A') = \A'\. We make an assumption that a move can only be done if it gives non-equivalent 
complex. We formulate the following lemma without proof. 

Lemma 14 If A 2 = /t = then the Markov chain on each A' is reversible with respect to the uniform 
measure. On the class of simplicial complexes this component coincides with the whole class. 

Proof. Reversibility is verified via the condition ir a \ a p = c(A') = 7I 7?'Va if ^(3a = X a /3 = X. 
The following example shows that large time and large N limits are not interchangeable, that is 

lim lim ^ lim lim 

t^co N— »oo AT— »oot— >oo 

for local quantities. This the simulation is slow and dangerous in this case. Consider the sequence of 
such chains ^ N \t) having the embedded state spaces 

... c A N c A N+1 c ... 

Take a vertex v at time and consider random variables qi N ^ (t) - number of edges at v at time t. We 
have L = V = and it could be natural to think that q v y — > 3. But the following argument 
shows more complicated situation. 

Lemma 15 Consider the class of simplicial complexes. As N — ► oo the limiting process exists and is 
the random walk on [3, oo) with transition rates Aj^+i = A^j-i = Xi. Thus P(qi N \t) = k) —>n->oo 0. 

Proof. For fixed iV the process q\j ^(t) is Markov with state space 3, ...,N with rates Aj^+i — 
Aj ; j_i = Xi. In fact, each edge incident to v can be changed to a transversal and, for each triangle 
containing v, its edge not containing v can be erased by GV-move, this will give one more incident 
edge. The limiting random walk is null recurrent and thus big fluctuations in it occur until it reaches 
equilibrium for fixed N. 

Similar proof does not hold for other classes of complexes. 

Now consider Markov chains where the only transitions are A- moves. To get ergodic chains we 
change the generator which produces jumps. Now the jumps are produced by any vertex i with rates 
A or /t. For fixed i with rate A take randomly (that is with probability q^ 1 ) one of the edges on the 
boundary of St(i) and do the A-move corresponding to this edge. Let \x be the rate of the inverse 
A-move at vertex i, also for each possible vertex v of degree 4 on dSt(i) with equal probability we 
take one pair of triangles (on the right hand side of the A-move) and do the inverse A-move. Once 
the vertex i appeared it can disappear afterwards. Let t(i) the time when vertex i appeared. 
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Theorem 13 If A > fj, then qi(t) — > oo mi/j positive probability. If X < fi then the vertex disappears 
a.s. and Eqi{t) is uniformly bounded. 

Proof. Let for each vertex v a(t) = a v (t) be number of vertices j on dSt(v) with qj = 4, let 
b(t) = Qv{t) — a(t). Fix vertex i. If v € dSt(i) then u + 1 is the next vertex on dSt(i) in the clockwise 
direction. 

Consider the process hit)) and for fixed configuration outside St(i) write down its infinites- 

imal jumps in Z?. For a direct and inverse A-move there can be only three possibilities: 




Figure 10: Proof of the theorem 



1. Two edges (marked 1 on the figure) appear on some link ((« — 1, v) on the figure). Thus here the 

transition is a, b — > a + 1, b with rate 2A- = 2A. Here and further factor 2 because the same move 

q 

can be produced also by the opposite vertex. Inverse move a, b — > a, b + 1 with rate 2//| = 2/j,; 

2. This move is produced by vertex v G dSt(i) (dotted edges 2 on the figure), the new vertex 
appears on the edge (i,v + 1). It produces a change in the vector (a, b) only if q v+ \ / 4. Thus 
here a,b — ► a + 1, 6 — 1 with rate 2Aa— . Inverse move gives the jump a, 6 — > a — 1, b with rate 

3. Next move is also produced by vertex v (edges 3 on the figure). q v can be transformed 4 — > 
3,5 ^ 4. 

In fact we do not need rates for 2 and 3: note only that these jumps conserve qi = a + b. Assume 
first \x > A. Then the embedded process f(n) = qi(t n ) = a(t n ) + b(t n ), where t n are the jump moments, 
satisfies the following inequality 

E(f(n + l)\f(n))<f(n)-e 

for some fixed e > 0. By the submartingale techniques (see, for example, ^8|) we have the proof. In 
the opposite case we have 

E{f(n + 1)\ f(n))> f(n) + e 
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and again the techniques of [48 works. 

It seems plausible that if A < [i then for all sequences t(i) — > oo the process qi{t) tends to some 
proper distribution if s = t — t(i) is fixed. If A <C /x it can be proved. On the contrary for the critical 
case fx = A random variables qi{t) fluctuates as for the Brownian motion. Compared with the results 
in the previous section this gives argument that we do not get the physical invariant measure here. 
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